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Abstract 

In this paper, we study biharmonic maps into Sol and Nil spaces, two model 
spaces of Thurston's 3-dimensional geometries. We characterize non-geodesic 
biharmonic curves in Sol space and prove that there exists no non-geodesic 
biharmonic helix in Sol space. We also show that a linear map from a Eu- 
clidean space into Sol or Nil space is biharmonic if and only if it is a harmonic 
map, and give a complete classification of such maps. 

1. INTRODUCTION 

In this paper, we work in the category of smooth objects, so manifolds, maps, 
vector fields, etc, are assumed to be smooth unless it is stated otherwise. 

A map if : (M, g) — > (N, h) between Riemannian manifolds is called a bihar- 
monic map if cp\fl is a critical point of the bienergy 

?2 I, „ n\_ 1 f l_ /,„M2. 



E*(<p,Sl) = - IrafoOrda: 

for every compact subset Q of M, where T2(<^) = Trace s Vd</? is the tension field 
of tp. Using the first variational formula (see [13]) one sees that ip is a biharmonic 
map if and only if its bitension field vanishes identically, i.e., 

(1) r 2 (ip) := - A^foO) - Trace^d^, r 2 (ip))dif = 0, 

where 

= -Trace 3 (V^) 2 = -Trace 9 (V^V v - V^„) 
is the Laplacian on sections of the pull-back bundle tp~ l TN and R N is the cur- 
vature operator of (N, h) defined by 

R N (X,Y)Z=[V^}Z-V^ Y] Z. 

Note that t 2 (</?) = —jvfofa)), where J v is the Jacobi operator which plays an 
important role in the study of harmonic maps. 
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Clearly, any harmonic map is biharmonic, so harmonic maps are a subclass of 
biharmonic maps. It is interesting to note that so far, apart from maps between 
Euclidean spaces (in which case any polynomial map of degree less than four is 
a biharmonic map), not many example of proper (meaning non-harmonic) bihar- 
monic maps between Riemannain manifolds have been found (see, e.g., [TO], [9], 
[TT] , and the bibliography of biharmonic maps [8] ) . 

Recently, some work has been done in the study of non-geodesic biharmonic 
curves in some model spaces. For example, [2] gives a complete classification 
of non-geodesic biharmonic curves in 3-sphere while in [3] it is proved that a 
non-geodesic biharmonic curve in Heisenberg group H 3 is a helix and the explicit 
parametrizations of such curves are given. For the study of biharmonic curves 
in Berger's spheres, in Minkowski 3-space, in Cartan-Vranceanu 3-dimensional 
space, and in contact and Sasakian manifolds see [T], [7], [I], and [B] respectively. 

In this paper, we first characterize non-geodesic biharmonic curves in Sol space 
(Theorem 12.31) and show that there exists no non-geodesic biharmonic helix in 
Sol space (Theorem 12. 71) . In the second part of the paper, we study linear bi- 
harmonic maps from Euclidean space R m into Sol space (R 3 ,^,,;) and Nil space 
(M. 3 ,gNu) using the linear structure of the underlying manifolds of Sol and Nil 
spaces. We show that such a linear map is biharmonic if and only if it is a har- 
monic map and give a complete description of such maps ( Theorems 13.31 and 13.41) . 



2. Biharmonic curves in Sol space 

Biharmonic curve equation in Frenet Frames. Let / C R be an open 
interval, and 7 : / — ► (N, h) be a curve, parametrized by arc length, on a 
Riemannian manifold. Putting T = 7', we can write the tension field of 7 as 
T (l) — Vy7', and the biharmonic map equation ([1]) reduces to 

V^T- R N (T, V T T)T = 0. 

A successful key to study the geometry of a curve is to use the Frenet frames 
along the curve which is recalled in the following 

Definition 2.1. (see, for example, \V2\ ) 

The Frenet frame {i 7 i}i=i ) 2,...,n associated to a curve 7 : / C R — > (N n ,h) 
parametrized by arc length is the orthonormalisation of the (n+l)-tuple 
{V ( J; ) d7(g)} fc=0 ,i,2, ..., n described by: 
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V^Fi = k x F 2 

vfi^ = -h^F^ + hF i+1 ,Vi = 2, 3, . . . , n - 1 

at 

where the functions {k\, k 2 , ■ ■ ■ , k n _i\ are called the curvatures 0/7 and V 7 is the 
connection on the pull-back bundle / -f~~ 1 (TN). Note that Fi = T = 7' is the unit 
tangent vector field along the curve. 

With respect to the Frenet frame, the biharmonic curve equation takes the 
following form. 

Lemma 2.2. (see, for example, Let 7 : / C K. — > (N n ,h)(n > 2) be a 
curve parametrized by arc length from an open interval of M into a Riemannian 
manifold (N n , h). Then 7 is biharmonic if and only if : 

hk[ = 

k'l — k\ — hk 2 2 + k 1 R(F ll F 2 , F x , F 2 ) = 

(2) { 2k' 1 k 2 + k 1 k' 2 + k 1 R{F 1 ,F 2 ,F 1 ,F 3 ) = 

k 1 k 2 h + hR(F 1 ,F 2 ,F u Fi) = 

I k 1 R(Fi,F 2 , Fi, Fj) = i = 5,...,n. 

Riemannian structure of Sol space. Sol space, one of Thurston's eight 
3- dimensional geometries, can be viewed as 1R 3 provided with Riemannian metric 
Qsoi = ds 2 = e 2z dx 2 + e~ 2z dy 2 + dz 2 , where (x, y, z) are the standard coordinates 
in M. 3 . Note that the Sol metric can also be written as: 

3 

ds 2 = y^o/ ® uj\ 
i=l 

where 

uj 1 = e z dx, oj 2 = e~ z dy, u 3 = dz, 

and the orthonormal basis dual to the 1-forms is 

_ z d d d 

e i = e tt, e 2 = e —, e 3 = — . 
ox ay oz 

With respect to this orthonormal basis, the Levi-Civita connection and the Lie 
brackets can be easily computed as: 

V ei ei = -e 3 , V ei e 2 = 0, V ei e 3 = e x 

(3) V e2 ei = 0, V e2 e 2 = e 3 , V e2 e 3 = -e 2 

V es ei = 0, V e3 e 2 = 0, V es e 3 = 0. 
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(4) [ei, e 2 ] = 0, [e 2 , e 3 ] = -e 2 , [ei, e 3 ] = e x . 

We adopt the following notation and sign convention for Riemannian curvature 
operator. 

(5) R(X, Y)Z = V X V Y Z - V Y V X Z - V [X , Y] Z, 
the Riemannian curvature tensor is given by 

(6) R(X, Y, Z, W) = g(R{Y, X)Z, W) = -g(R(X, Y)Z, W). 
Moreover we put 

(7) R a bc = R(e a , e b )e c , R a bcd = R(e a , e&, e c , e d ), 
where the indices a, b, c, d take the values 1, 2, 3. 

A direct computation using (jBJ), (BJ, (IS), (EJ), and (J7J) gives the following non- 
zero components of Riemannian curvature of Sol space with respect to the or- 
thonomal basis {ei, e 2 , 63} (we do not list those that can be obtained by symmetric 
properties of curvature): 

R121 = ~ e 2 , i?i3i = e3, i?i 22 = e±, 
R232 = 63, -R133 = — ei, i? 2 33 = — e 2 , 

and 

#1212 = -#(-R(ei,e 2 )ei,e 2 ) = -g(-e 2 ,e 2 ) = 1, 

(8) ^1313 = -g{R{ei, e 3 )ei, e 3 ) = -#(e 3 , e 3 ) = -1, 
-R2323 = -fi'( J R(e 2 ,e 3 )e 2 ,e 3 ) = -g(e 3 ,e 3 ) = -1. 

Biharmonic curves in Sol space. Let 7 : / — ► (M. 3 ,g,soi) be a curve on 
Sol space parametrized by arc length. Let {T, N, B} be the Frenet frame fields 
tangent to Sol space along 7 defined as follows: T is the unit vector field 7' 
tangent to 7, N is the unit vector field in the direction of VtT (normal to 7), 
and B is chosen so that {T, N, B} is a positively oriented orthonormal basis. 
Then, by Definition (12.11) . we have the following Frenet formulas 

( V T T = kN 

(9) < V T iV = —kT + tB 

{ S/ T B = -tN, 

where k = |VtT| is the geodesic curvature and r the geodesic torsion of 7. With 
respect to the orthonormal basis {ei,e 2 ,e 3 } we can write T = T x t\ + T 2 e 2 + 
T 3 e 3 , N = Nie ± + N 2 e 2 + N 3 e 3 , B = T x N = B x e x + B 2 e 2 + B 3 e 3 , and we have 
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Theorem 2.3. Let 7 : / — > (M?,gs i) be a curve parametrized by arc length. 
Then 7 is a non- geodesic biharmonic curve if and only if 



k = constant 7^ 
(10) { k 2 + r 2 = 2B 2 -l 

t' = 2N 3 B 3 . 



Proof. By §2§ of lemma (\2.2\i we see that 7 is a biharmonic curve if and only if 



kk! = 

k" - k 3 - kr 2 + kR(T, N, T, N) = 
2k' t + kr' + kR(T, N, T, B) = 0, 



which is equivalent to 



k = constant 7^ 

(11) { k 2 + t 2 = R(T, N, T, N) 

t' = —R(T, N, T, B) 



since k 7^ by the assumption that 7 is non-geodesic. A direct computation 
using ([H]) yields 



R(T, N, T, N) = Elj,i, P =i TiNpTiNjRipij 

= TxNzTiNiRnn + T 1 N 2 T 2 N 1 R X22l + T 2 N 1 T 2 N 1 R 2X2l + T 2 N 1 T 1 N 2 R 2112 

+T 1 N 3 T 1 N 3 R 1313 + T 1 N 3 T 3 N 1 R 1331 + T-^T^R^ + T 3 NxTiN 3 R 3113 

+T 2 N 3 T 2 N 3 R 2323 + T 2 N 3 T 3 N 2 R 2332 + T 3 N 2 T 3 N 2 R 3232 + T 3 N 2 T 2 N 3 R 3223 

= T 2 N 2 - TxT 2 N x N 2 + T 2 N 2 - T X T 2 N X N 2 

-T 2 N 2 + TxTgiVxiVa - T 2 N 2 + T X T 3 N X N 3 

-T 2 N 2 + T 2 T 3 N 2 N 3 - T 2 N 2 + T 2 T 3 N 2 N 3 

= [T X N 2 - T 2 N X ) 2 - (T 3 Nt - T X N 3 ) 2 - (T 3 N 2 - T 2 N 3 ) 2 

— ts 3 - tt x - tt 2 — zr> 3 - 1 

(by T x N = B ,T x B = -N, B 2 + B 2 + B 2 = 1), 
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and 

R(T, N, T, B) = E- iAp=1 TiNpTiBjRipij 

= T1N2T1B2R1212 + T1N2T2B1R1221 

+T2N1T2B1R2121 + T2N1T1B2R2112 

+Ti N3T1B3R1313 + ^Ns^BiRi^i 

+T 3 7V 1 T 3J B lJ R3i3i + T 3 N X T X B 3 R 3113 

+T2N3T2B3R2323 + T2N3T3B2R2332 

+T3N2T3B2R3232 + T3N2T2B3R3223 

= T?N 2 B 2 - T 1 T 2 N 1 B 2 + N X B X - T 1 T 2 N 2 B 1 

-T 2 N 3 B 3 + T X T Z N X B 3 - T^N 1 B 1 + T 1 T 3 B 1 N 3 

-T^N 3 B 3 + T 2 T 3 N 2 B 3 - T^N 2 B 2 + T 2 T 3 B 2 N 3 

= (T 1 B2-T2B 1 )(T 1 N 2 -T 2 N l ) 

—(T 3 Nx — T 1 N 3 )(T 3 B 1 — TxB 3 ) 

-(T 3 N2-T2N 3 )(T 3 B 2 -T 2 B 3 ) 

= -N 3 B 3 + N 2 B 2 + N 1 B 1 = -2N 3 B 3 

(by T x N = B ,T x B = —N, + + iViS 3 = 0), 

these, together with Equation (II ip . complete the proof of the theorem. □ 

As an immediate consequence we have 

Corollary 2.4. Let 7 : I — ► (M?,gs i) be a non- geodesic curve parametrized by 
arc length. If B 3 = 0, then 7 is not biharmonic. 

Corollary 2.5. Let 7 : I — ► (M?,gs i) be a non-geodesic curve parametrized by 
arc length. If B 3 is constant and N 3 B 3 7^ 0, then 7 is not biharmonic. 

Similar to the terminology used for curves in M 3 , we keep the name helix for a 
curve in a Riemannian 3-manifold having constant both geodesic curvature and 
geodesic torsion. With this terminology, we can use Equation (fill to deduce the 
following 

Corollary 2.6. Let 7 : / — > (M. 3 ,gs i) be a non-geodesic biharmonic helix 
parametrized by arc length, then 

{B 3 = constant 7^ 
N 3 = 
k 2 + r 2 = 2B 3 - 1. 

Non-geodesic biharmonic helices in 3-dimensional sphere S 3 , in Heisenberg 
group H 3 , and in Cartan-Vranceanu 3-dimensional space have been studied in 
[2], [3], and [1] respectively. In contrast to the situations in those 3-dimensional 
spaces where there are rich examples of such curves our next theorem shows that 
there exists no such curves in Sol space. 
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Theorem 2.7. There exists no non- geodesic biharmonic helix in Sol space. 

Proof. Suppose that 7 : I — > (M. 3 ,g So i) is a non-geodesic biharmonic helix 
parametrized by arc length. We shall derive a contradiction by showing that 
7 must be a geodesic. We can use to compute the covariant derivatives of the 
vector fields T, N, B as: 

( V T T = (T[ + 7\T 3 )ei + (T 2 ' - T 2 T 3 )e 2 + (T| - T 2 + T£e 3 

(13) I V T N = (N[ + HJV 3 )ei + (A^ - T 2 A^ 3 )e 2 + {T 2 N 2 - + N£e 3 
{ W T B = (B[ + T 1 B 3 )e 1 + (B' 2 - T 2 B 3 )e 2 + (T 2 B 2 - T X B X + B' 3 )e 3 . 

It follows that the third components of these vectors are given by 

r (V T T,e 3 >=T 2 2 -T 2 + T 3 ' 

(14) I (V T N,e 3 ) = T a iV 2 -T 1 ;V 1 + JV£ 

( (V T B,e 3 ) = T 2 B 2 -T l B 1 + B' 3 . 

On the other hand, using Frenet formulas (Q we have, 

( (V T T,e 3 ) = kN 3 

(15) I (V T N,e 3 ) = -kT 3 + rB 3 

( (V T B,e 3 } = -rN 3 . 

Since 7 is assumed to be a non-geodesic biharmonic helix, we have, by Corollary 
N 3 = 0, B 3 = constant. These, together with Equations ffT4"|) and ffT5"]) . give 



Tf — Tf + T3 = 

( 1 (i i <( T 2 A^ 2 - TxTVi = -A;T 3 + tB 3 

T 2 B 2 - T 1 B 1 = 0. 

Noting that T x B = —N we also have 

(17) T 2 B X - T X B 2 = N 3 = 0. 

Thus, we have 

T 2 N 2 - TliVi = -fcT 3 + r5 3 , (1) 
2f-2? + I^ = 0, (2) 
T 2 B 2 -T 1 B 1 =0, (3) 
T 2 B l -T l B 2 = $. (4) 

Case A: 7^ T 2 2 . In this case, Equations (3) and (4) in System ffTg]) viewed as 
equations in Bi and B 2 has a unique solution B\ = B 2 = 0. This implies that 
T 3 = (N x S, e 3 ) = 0. Substitute this into (2) of System JTHJ) we have T 2 = T|, 
a contradiction. Thus, we must have 

Case B: T 2 = T 2 2 . In this case, Equation (2) of System ffTgj) implies that T 3 = 
constant. To understand the meaning of this, we represent the unit tangent vector 



(18) 



8 YE-LIN OU AND ZE-PING WANG 

T as T = sinacos/?ei + sinasin/3e 2 + cosae 3 , where a = a(s), (5 = f3(s). 
With this representation, T 3 = constant implies that cos a = constant and hence 
a(s) = «o, a constant. This, together with T 2 = T|, gives 

(19) sin a (cos/5 ± sin (3) = 0. 

If sincto = 0, then we have T\ = T 2 = 0, and it follows from the first equation of 
( ITBl that VtT = which means that 7 is a geodesic, a contradiction. Thus, we 
must have sina 7^ 0, which, together with ( JT91 . implies that 

cos/? = ±sin/3 = 

and hence, 

\/2 

(20) ri = ±r 2 = ±-ysiiiao. 
We use the first equation of ( {TBI again to get 

VtT = sinoiocosa!o(±-z-ei =F -7^2) = kN 

which yields 

(21) N 1 = ±y^,N 2 = T y^- 

since = IVtT 1 ! = | sin ceo cos ao I • By the assumption that 7 is non-geodesic, we 
may assume, without loss of generality, that sinaocosao > 0, so 

(22) k = sin a cos ct . 

Using Equations ( l20i) . (12ip and the fact that B = T x N we have 

2?i = constant, 

(23) £> 2 = constant, 

B 3 = nN 2 - T 2 N X = =f sin cjq. 

It follows from (3) of (TT5|). and the third equation of (Q3]) that 

(24) r 2 = |V T 5| 2 = (T 2 + Tl)B\ = sin 4 a . 
Substituting (1221. (123)1 and (j22j) into the third equation in ([JJ) we have 



which implies 
and hence 



sin 2 a cos 2 a + sin a = 2 sin 2 a — 1 , 



sin 2 ckq = 1, 



cos 2 a 
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It follows that cos«o — from which and (1221) we conclude that k = 0, i.e., 7 is 
a geodesic, a contradiction. 

Combining the the results in Cases A and B we complete the proof of the 
theorem. □ 

3. Linear biharmonic maps into Sol and Nil spaces 

In this section, we first derive the biharmonic map equation in local coordinates 
and then we use it to classify linear biharmonic maps from Euclidean space into 
Sol and Nil spaces. In the rest of the paper, we adopt the Einstein summation 
convention that a repeated upper and lower index means the summation of that 
index over its range that is understood from the context. 
3.1 Biharmonic map equation in local coordinates. 



Lemma 3.1. Letip : (M m ,g) — > (N n ,h) with (p(x\ . . . , x m ) = (^(x), . . . , <p n (x)) 
be a map between Riemannian manifolds. With respect to local coordinates (x l ) in 
M and (y a ) in N, if is biharmonic if and only if it is a solution of the following 
system of PDE's 

(25) g^r* + r>ff^ + ^(r>?fy + r^^K P 

-I* (t£ + - r»^ a J = 0, a = 1, 2, . . . , n. 

Proof. Let gij and (resp. h a p and f ^) denote the components of metric and 
the connection coefficients in the domain (resp. the target) manifold with respect 
to the chosen local coordinates and the natural frame {^1} (resp. {^}). The 
bitension field of ip can be computed as 



r%) = <p V^VV - V£ _al (r(<p)) - Trace^c^r^))^ 
= 9 ij V%V% (r(cp)) - ^r* V^ e (rfo)) 

dx^ dx3 dx^ 

(26) -!?R{M^),«<p))M^) 



d d 

dx i ax k oy uy 1 j 

A direct computation gives 



g" (t( V )) - T%V% (rl») - ^r { " t{v)) . 



(27) V A <r<„)) = V^r" A) = (r; + T ^ A ) A 
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V^V^Jr^)) = Vn L VV(^7^ 7 )=V^frf + rVff^ ° 

dx % dx3 dx 1 dx3 



8y"> W -&V* ^ A *7 8y* 



(28) 
and 



/ - 8 - \ 8 



(29) R{i-M^))i- = T»W 



8 8-8 
dr ,T{ip)) W = rURr7pau dy~° 



Substitute Equations (l27jl . (1251). and (123]) into (EE} we have 



(30) -4(r fc CT + T>fry - rVi^oJ^, 



from which the lemma follows. □ 

When the domain manifold is a Euclidean space then we have 

Corollary 3.2. Let <p : M m — ► (iV n , ft) wt/i y?(x\ . . . , x m ) = (^(x), f n {x)) 
be a map from a Euclidean space into a Riemannian manifold. Then ip is bihar- 
monic if and only if it is a solution of the following system of PDE's 

(31) Ar CT + (Vr°, V^)f^ + (V/, V(r Q f ^ )) 



a/3 "T- \ v ^ , v v J- aj g> 

+ (V/, Vip»)T a t v ap t % - r"(V<^, V/)^ a , = 0, a = 1, 2, . . . , n. 



Linear biharmonic maps into Sol space. Let (M 3 ,^,,;) denote Sol space, 
where the metric can be written as gsoi = e 2y (dy 1 ) 2 + e~ 2y ' (dy 2 ) 2 + (dy 3 ) 2 with 
respect to the standard coordinates (y 1 , y 2 , y 3 ) in M 3 . Then a direct computation 
gives the following components of Sol metric and the coefficients of the connection: 



9n = e 2y , g 2 2 = e 2y , g 33 = 1, all other g id = 0; 

e 2 y\ g 33 = 1, all other g ij = 0; 



,11 _ ^-2y 3 „22 _ 2j, 3 33 
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(32) 



r 1 

11 




0. 


r 3 - - 
ii — 


e 2yz 


f 1 

12 


= o, r 2 2 




o,r 3 2 = 


0; 


f 1 

13 




— 


o,f 3 3 = 


0; 


f 1 

21 


= 0,^ 


— 


o,rii = 


0; 


f 1 


- r 2 - 

ZZ 


0. 


ZZ 


-2y :i . 


f 1 

1 23 


= o,rf 3 




-1 f 3 
- 1 ! 1 23 


= 


1 31 


= i,r§! 




n r 3 - 

u > 1 31 - 


0; 


1 32 


= o,r§ 2 




"l,f 3 2 


= 


1 33 


= o,f| 3 




o,r| 3 = 


0. 



By our convention of curvature operator and the following notation for the com- 
ponents of the Riemannian curvature 



(33) 

we have 
(34) 



d d , d 



R 



k i] 



_d_ 

dy l 



R 



k i.i 



d d 

j^l | -p/ -p£ -p/ -pt 

Qyi kj ~ Qyj 1 ki i" 1 it 1 fej — 1 ji 1 ki 



A straightforward computation using (J32l) and (l34l gives the following compo- 
nents of the Riemannian curvature of Sol space: 



pi 

-"■221 



-2r 



(35) 



p3 
-"-131 



p2 
rl 121 

-e 2 ^,i? 3 



) ^3 

e 2 y\R 2 



331 
j 

112 : 



: 1) -^212 ~~ e 
e > rL 332 " 



2y 3 Dl 
> -"-313 

- 1 P 2 
" rL 323 



-i; 

-1; 



232 



-2y" £ 



113 



o2y 3 d3 
> -^223 



Now we are ready to prove the following classification theorem for linear bihar- 
monic maps into Sol space. 



Theorem 3.3. Let ip:R m — ► (M 3 , g So i) with 

( x 1 \ 

x 2 

\ X m j 

i.e., <p(X~) = (A 1 X t , A 2 X l , A 3 X f ) be a linear map into Sol space, where A % denotes 
the row vectors of the representation matrix. Then, ip is a biharmonic map if and 
only if it is a harmonic map, which is equivalent to either (i) A 3 = 0, \A Y \ 2 = 
\A 2 \ 2 , or, (ii) A 3 ^ 0, A 1 = A 2 = 0. 



ip(x) 



a\ 


«2 • ■ 


■ a 1 


a 2 


a 2 ■■ 


■ °?m 


a 3 


4 ■■ 


■ <& 
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Proof. With respect to the standard Cartesian coordinates (x l ) in IR m and (y a ) 
in IR 3 , the tension field of ip is given by 



r(ip) : = Trace 9 (Vd^) e T^TN) 

_ _<T d 

■=i 

cr 



J=l 



where and in the sequel, A a ■ A 13 denotes the inner product and \A a \ the norm of 
the vectors in Euclidean space. 

Putting r((p) = t° and substituting ( 1321) into Equation (1361) we find the fol- 
lowing components of the tension field of ip 



r 1 = 2A 1 -A 3 , 
(37) r 2 = -2A 2 ■ A 3 , 

r 3 = \A 2 \ 2 e- 2 y 3 - \A l \ 2 e 2 y 3 , 



where and in the sequel y 3 = A 3 X f . 
A further computation gives 



(A^+(Vr«,V/)f^) A 
(At* + ^ • Vr<T A = 



(38) 


-2A 1 


■A 3 (|A 2 | 2 e- 2y3 + 


(A 1 ) 




+2A 2 


•A 3 (|vl 2 | 2 e- 2y3 + 


1^1 



,,3, 



<9 



dy 2 



+4|A 3 | 2 (|A 2 | 2 e- 2 ^ 3 - lA 1 ^) A 
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(39) 



d_ 



= (^-Vr«r^ + r^-vr^)A 

= -2A 1 ■ A 3 (\A 2 \ 2 e~ 2y3 + \A l \ 2 e 
+2A 2 ■ A 3 (\A 2 \ 2 e- 2y3 + \A 1 rr 



dy 1 



} dy 2 



+ [-A(A 1 ■ A 3 ) 2 e 2y + A{A 2 ■ A 3 )~e 



2^-2^ 



l dy 3 



(40) 



and 



= [2A 1 ■ A 3 \A 3 \ 2 + A 1 ■ A 3 \A 2 \ 2 e- 2yi - 3A 1 ■ A 3 \A 1 \ 2 e 2y3 
-2A 2 ■ A 3 A 1 ■ A 2 e- 2y3 }^ + [-2A 2 ■ A 3 \A 3 \ 2 - A 2 ■ A 3 \A l \ 2 e 2y3 
+3A 2 ■ A 3 \A 2 \ 2 e~ 2y3 + 2A 1 ■ A 3 A 1 ■ A 2 e 2y3 }^ + \-2{A x ■ A 3 fe 2y3 



+|A 1 | 4 e 42/ + 2(A 2 ■ A 3 ) 



\A 



9y 

2|4„-42/ 3 l_a 



(V^, V^)R^ au ^ = AP ■ A^Rl^ 



[-3A 1 • A 3 \A 2 \ 2 e~ 2yA - 2A 2 ■ A 3 A 1 ■ A 2 



-2y A 



(41) +A 1 ■ A 3 \A 1 \ 2 e 2y + 2A 1 ■ A 3 \A 3 \ 2 ]^ + [2 A 1 ■ A 3 A 1 ■ A 2 e 
+3 A 2 ■ A 3 \A 1 \ 2 e 2y3 - A 2 ■ A 3 \A 2 \ 2 e~ 2y3 - 2 A 2 ■ A 3 \A 3 \ 2 }^ 
+ [-|A 1 | 4 e 4y3 + \A 2 \ A e~ iy3 - 2(A 1 ■ A 3 ) 2 e 2y3 + 2(A 2 ■ A 3 ) 2 e~ 2y3 )}-7^ 



dy 3 ' 



It follows from Equations (j5gj|. (1391. (HOI). (TO and Corollary E2 that the linear 
map (p is a biharmonic map if and only if 



(42) 



-ZA 1 ■ A 3 \A x \ 2 e 2y3 =0 
8A 2 ■ A 3 \A 2 \ 2 e~ 2y3 =0 
4(|yl 2 | 2 |A 3 | 2 + (A 2 ■ A 3 ) 2 )e~ 2y3 - 4(|A 1 | 2 |A 3 | 2 + (A 1 ■ A 3 ) 2 )e 2y " 



+2\A 



i|4 p v _ 2\a 2 \ a p-^ 



Solving System of equations (|42[) we have either (i) A 3 = 0, l^ 1 ) 



: 0. 



or 

(ii) A 3 ^ 0, A 1 = A 2 = 0. It follows from Equation (l3Tj) that in both cases the 
tension field vanishes identically, i.e., <p is also harmonic. Therefore, we obtain 
the theorem. □ 



Linear biharmonic maps into Nil space. Let (R , g^u) denote Nil space, 
where the metric with respect to the standard coordinates {y l ,y 2 ,y 3 ) in M 3 can 
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be written as g^u = (Ay 1 ) 2 + (dy 2 ) 2 + (ch/ 3 — y 1 dy 2 ) 2 . Then an easy computation 
gives the following components of Nil metric and the coefficients of its connection: 

9X1 = 1, 012 = 013 = 0, 922 = 1 + G/) 2 , #23 = -V 1 , 033 = 1| 

11 = 1, 12 = 13 = 0, g 22 = 1, s 23 = y 1 , ^ 33 = 1 + (y 1 ) 2 ; 



(43) 



r 1 
1 ii 


= f 


2 - n 
ii — u i 


r 3 - n- 
1 ii — u i 




f 1 

1 12 


= 


f 2 - 

1 12 — 


y 1 f3 . 

2 ' 1 12 ~ 


. (y 1 ) 2 -! 

2 


1 13 


= 


f 2 - 
1 13 — 


1 p3 
2' 1 13 


_ V 1 . 
2 ' 


1 21 


= 


f 2 - 
1 21 ~~ 


y 1 f3 _ 

2 ' 1 21 ~ 


2 


1 22 




V l S 2 22 


= o,r 3 2 


= 0; 


1 23 


1 

2 


V 2 - 

■> 1 23 ~~ 


o, r| 3 = 


0; 


f 1 
1 31 


= 


f 2 - 
1 31 — 


_i f 3 

2' 1 31 


_ y 1 . 

2 ' 


1 32 


1 

2 


r 2 - 

i 1 32 — 


o, f 3 2 = 


0; 


1 33 


= 


r 2 - 

1 33 — 


o,r 3 3 = 


0. 



A further computation using ( 134")) and ( 143i) gives the following components of 
the Riemannian curvature of Nil space: 

fa 1 ) 2 51 _ _yi 51 _ 3 fa 1 ) 2 
5 ^213 ~~ 4 ' lx 221 ~ 4 4 ' 

4 ' ^313 



(44) 



pi 

rt 212 



- 5 + 
4 T 



'213 
pi _ 



pi 
rl 312 



1 51 _^ Dl 
4 ) rt 321 ~~ 4 > -"-3 



i? 2 



112 — 
-"-232 



4 ' * "321 
^121 



3 p2 



,1 



323 



p3 
- rl 223 



R3 

-"-112 " 



y\R 3 



, .ft 



3 

232 



113 — 4> -"-121 

R3 _ 

-"-323 " 



331 

_3 p2 _ 
4 ) rl 223 " 
_ 1 52 
4 ' -"-332 



-y\R 3 



4 



131 



4 ' -"-332 



r_ 

4 ' 
_ 1. 

4' 

ai 

4 ' 
_ 1. 

4' 
_ 1 

4' 

4 ' 



Theorem 3.4. Let ip : 



(p(x) 



\gmi) with 



a\ 


«2 • ■ 


■ a 1 


a 2 


«| • ■ 


■ a 2 


a 3 


a% ■ ■ 


■ a 3 



( x 1 \ 

x 2 



\x m J 

i.e., tp(X) = (A 1 X t , A 2 X l , A 3 X t ) be a linear map into Nil space, where A 1 denotes 
the row vectors of the representation matrix. Then, ip is a biharmonic map if and 
only if it is a harmonic map, which is equivalent to either (i) A 1 = 0, A 2 ■ A 3 = 0, 
or, (ii) A 1 = 0, A 2 = 0, or, (Hi) A 2 = 0, A 1 ■ A 3 = 0. 



Proof. Taking the standard Cartesian coordinates (x % ) in IR m , (y a ) in 
substituting (JI3"]) into flHol) we find the tension field of p> to be 



and 
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dy" 

+[A 1 -A 2 y 1 -A 1 -A 3 ]- S 



J dy 2 

+ [A^A\y^-A'.A^-A^A*]JL 



or, 



r 1 = -|A 2 |V + A 2 ■ A 3 , 
(46) t 2 = A 1 ■ A 2 y l - A 1 ■ A 3 , 

t 3 = A 1 - A 2 {y 1 ) 2 - A 1 ■ A 3 y l - A 1 ■ A 2 , 

where and in the sequel y 1 = A 1 X t . 
A straightforward computation yields 



(47) Ar^ = 2A^A 2 \A^\ 2 ^, 



(48 



(W, V</>f = A? • Vt°F 

= [-IA 1 ■ A 2 (\A'\ 2 + |A 2 | V + \A l ■ A 3 (\A l \ 2 + \A 2 \ 2 )]^ 2 
+ • AWA^ 2 + \A 2 \ 2 )(y 1 ) 2 + \& ■ A 3 (\A^\ 2 + |A 2 | V 

+^.A 2 (\A 2 \ 2 -\A^\ 2 )]^, 



(V^,V(r-f^))^ 



(49) 



a/3 J/ dy" 

[-(A 1 ■ A 2 )^ 1 + (A 1 ■ A 2 )(A* ■ A 3 )}^ 



(A 13 ■ Vr a r^ + T a A^ ■ vr CT ^ 9 



+ [-A l ■ A 2 \A 2 \ 2 y l + \(A l ■ A 2 )(A 2 ■ A 3 ) + \A l ■ A 3 \A 2 \ 2 }^ 
+ [-\A l ■ A 2 \A 2 \ 2 (y 1 ) 2 + (A 1 ■ A 2 ){A 2 ■ A 3 )y l 
+A i . A 3 \A 2 \ 2 y l + \A X ■ A 2 \A 2 \ 2 - \{A X ■ A 3 )(A 2 ■ A 3 )}^, 
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T Si Si af3 L vpdy <j 



21 v 



+\{A 2 ■ A 3 )V - \{A l ■ A 2 fy l + \\A 2 \ 2 \A 3 \ V + j\A [ 
+\(A l ■ A 2 )(A l ■ A 3 ) - \A 2 ■ A 3 \A 3 \ 2 - \A 2 ■ A 3 \A 2 \ 2 ]^ 
+ [-\A l ■ A 2 \A 2 \ 2 (y 1 ) 3 + \A l ■ A 3 \A 2 \ 2 (y 1 ) 2 



12,2 1 



(50) +i(A 2 • A 3 )(A 1 ■ A 2 )^ 1 ) 2 - \A X ■ A 2 \A 
+ \A l ■ A 2 \A l \ V - iA 1 • A 2 |A 3 |V - A (A 1 • A 3 )(A 2 • A 3 )y l 
+\{A l ■ A 2 )(A 2 ■ A 3 ) + \A l ■ A 3 \A 3 \ 2 - \A l ■ A 3 ^ 1 ] 2 ]^ 

+ [-\A l ■ A 2 |A 2 | V) 4 + i^ 1 • A*\A 2 \\y x Y + \{A 2 ■ A 3 ){A l ■ A 2 )(y 1 ) 3 
+ \A X ■ A 2 \A 1 \ 2 (y 1 ) 2 - \(A 2 ■ A 3 )^ 1 ■ A 3 )^ 1 ) 2 - ^A 1 ■ A 2 \A 3 \ 2 (y 1 ) 2 
-IA 1 ■ A 3 |A 2 |V + lA 1 ■ A 3 |A 3 |y - \A l ■ A 3 \A l \ 2 y 1 
+\{A 2 ■ A 3 ){A l ■ A 3 ) + \A l ■ A 2 \A 2 \ 2 - \A l ■ A 2 ^ 1 ] 2 ]^ , 

and 

r"<V^, VipP)^^ = A? ■ A"T"R% aul ^ 

= [A|A 2 | V) 3 - \A 2 ■ A 3 \A 2 \ 2 (y 1 ) 2 - \{A X ■ A 2 ) 2 y x + \{A 2 ■ A 3 )V 
-l\A 2 \Y + i|A 2 | 2 |A 3 |y + ^A 1 ■ A 2 )(A l ■ A 3 ) + \A 2 ■ A 3 \A 2 \ 2 
-\A 2 ■ A 3 \A 3 \ 2 }^ + \-\A^ ■ A 2 \A 2 \ 2 ( y y + \A X ■ A 3 \A 2 \ 2 {y 1 ) 2 
+ \(A 2 ■ A 3 )(A l ■ A 2 )(y 1 ) 2 + \A X ■ A 2 \A l \ 2 y l + A 1 ■ A 2 \A 2 \ 2 y l 

(51) -\{A 2 ■ A 3 ){A l ■ A 3 )y x - \A l ■ A 2 \A 3 \ 2 y l - f-A 1 • A 3 \A 1 
-(A 1 ■ A 2 )(A 2 ■ A 3 ) + \A X ■ A 3 \A 3 \ 2 }^ + \-\A x • A 2 \A 2 \ 2 (y 1 ) 
+ IA 1 ■ A 3 \A 2 \ 2 (y 1 ) 3 + \(A 2 ■ A 3 ){A l ■ A 2 ){y 1 ) 3 + \A l ■ A 2 \A 1 \ 2 (y 1 ) 
+A 1 . A 2| A 2|2( yl )2 _ l (A 2 . A 3)( A i . ^(y 1 ) 2 - \A l ■ A 2 \A 3 \ 2 (y 1 ) 2 

-fA 1 • A 3 \A l \ 2 y l - (A 1 ■ A 2 )(A 2 ■ A^y 1 + \A l ■ A 3 |A 3 |y 
+ \A^ - A 2 \A'\ 2 + \A^ - A 2 \A 2 \ 2 }^. 

By Equations QM, Pty. (l5U]h (15II) and Corollary EJ we conclude that ip 
is biharmonic if and only if 



-(A X .A 2 )V + \A 2 \Y 
1 } +{A l ■ A 2 ){A l ■ A 3 ) - A 2 ■ A 3 \A 2 \ 2 = 0, 



-A 1 ■ A^A^y 1 - 3A l ■ A 2 \A 2 \ 2 y l 
(53) +A 1 ■ A 3 \A 2 \ 2 + A 1 ■ A^A 1 ] 2 

+2(A X ■ A 2 )(A 2 ■ A 3 ) = 0, 



2 

J - 1 

2 
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and 

-A 1 ■ A 2 \A 1 \ 2 (y 1 ) 2 - 3A 1 ■ A 2 \A 2 \ 2 (y 1 ) 2 + A 1 ■ A 3 \A 1 \ 2 y 1 
(54) +A 1 ■ A 3 \A 2 \ V + 2(A 1 ■ A 2 )(A 2 ■ A 3 )y Y 

+A 1 ■ A^A 1 ] 2 + A 1 ■ A 2 \A 2 \ 2 = 0. 

To solve the System of biharmonic map equations (1521) . ( 1531) . and (1541) we consider 
the following two cases. 

Case I: A 1 = 0. Noting that y 1 = A 1 X t = we use Equation ( J52l) to have 

-A 2 -A 3 \A 2 \ 2 = 0, 

which implies either A 2 ■ A 3 = 0, or \A 2 \ 2 = i.e. A 2 = 0. 

Substituting A 1 = 0, A 2 ■ A 3 = or A 1 = 0, A 2 = into pS) and {54} we see 
that they are both solutions of biharmonic map equations. 

Case II: A 1 ^ 0. Note that in this case, y l = A^-X 1 ^ 0. We can view sys- 
tem of biharmonic map equations ( 1521) , f l53i) , and f l54l) as a system of polynomial 
equations iny . It is not difficult to check that A 2 = 0, A 1 • A 3 = is the only 
solution in this case. 

Combining Case I and II we obtain the last statement of the theorem. A direct 
checking using the tension field equation (T4"5"]) we see that all these maps are also 
harmonic maps. Thus, we complete the proof of the theorem. □ 
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